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A finite rank separable approximation for the quasiparticle RPA calculations with 
Skyrme interactions that was proposed in our previous work is extended to take into 
account the coupling between one- and two-phonon terms in the wave functions of ex- 
cited states. It is shown that characteristics calculated within the suggested approach are 
"q ■ in a good agreement with available experimental data. 

1. INTRODUCTION 

Many properties of the collective nuclear excitations can be described within the random 
phase approximation (RPA) jll'2l3l4j . The most consistent models employ the Gogny's 
j3] or Skyrme-type [U] effective interactions which can describe the ground states in the 
framework of the Hartree-Fock (HF) approximation and the excited states within the 
RPA. Such models are quite successful for predicting nuclear states properties |7|8|9|1U|TT] . 

Due to the anharmonicity of vibrations there is a coupling between one-phonon and 
more complex states |2|4j and the complexity of calculations beyond standard RPA in- 
creases rapidly with the size of the configuration space, so one has to work within limited 
spaces. Using simple separable forces one can perform calculations of nuclear charac- 
teristics in very large configuration spaces since there is no need to diagonalize matrices 
whose dimensions grow with the size of configuration space. For example, the well-known 
quasiparticle-phonon model (QPM) j3] can do very detailed predictions for nuclei away 
from closed shells |T2l. 

The possibility for such a simplification was the motivation for proposing in our previous 
work a finite rank approximation for the particle-hole (p-h) interaction resulting from 
Skyrme-type forces. Thus, the self-consistent mean field can be calculated in the standard 
way with the original Skyrme interaction whereas the RPA solutions would be obtained 
with the finite rank approximation to the p-h matrix elements. It was found that the 
finite rank approximation reproduces reasonably well the dipole and quadrupole strength 
distributions in Ar isotopes |13j. 
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Recently, we extended the finite rank approximation for p-h interactions of Skyrme type 
to take into account pairing ^3]. In this paper we generalize our approach to take into 
account a coupling between the one- and two-phonon components of wave functions. As 
an application we present results of our first calculations for the quadrupole and octupole 
states in 112 Sn. 



2. METHOD OF CALCULATIONS 

We start from the effective Skyrme interaction [H] and use the notation of Ref.[JJ)] con- 
taining explicit density dependence and all spin-exchange terms. The single-particle spec- 
trum is calculated within the HF method. The continuous part of the single-particle 
spectrum is discretized by diagonalizing the HF hamiltonian on the harmonic oscillator 
basis 16]. The p-h residual interaction V res corresponding to the Skyrme force and includ- 
ing both direct and exchange terms can be obtained as the second derivative of the energy 
density functional with respect to the densitv[T7|. Following our previous paper [TB] we 
simplify V res by approximating it by its Landau-Migdal form. For Skyrme interactions all 
Landau parameters Fi,Gi,F t , G h with / > 1 are zero. Here, we keep only the I = terms 
in V Tes and in the coordinate representation one can write it in the following form: 

K es (ri, r 2 ) = Nq 1 [F (n) + Ga(n)ai(j 2 + (F^n) + (^(n)^)^] 5(n - r 2 ) (1) 

where Oi and Tj are the spin and isospin operators, and N = 2k F m* / 'ix 2 Ti 2 with k F and 
m* standing for the Fermi momentum and nucleon effective mass. The expressions for 
F , Go, F , G Q in terms of the Skyrme force parameters can be found in Ref . [To] . Because 
of the density dependence of the interaction the Landau parameters of Eq. (JJJ are functions 
of the coordinate r. 

The p-h residual interaction can be presented as a sum of N separable terms. Let us 
explain this procedure for making the finite rank approximation by examining only the 
contribution of the term Fq. In what follows we use the second quantized representation 
and V rP e can be written as: 



V r es = V 123 4 ■ afa%a A a 3 : (2) 

^ 1234 

where af (ai) is the particle creation (annihilation) operator and 1 denotes the quantum 
numbers (ni/ijimi), 



^1234 = J 0i(ri)02(r 2 )Ke S (ri,r 2 )03(ri)04(r2)dridr2, (3) 



V1234 Y^-I 2 (/i V/ ./3K/2 V, ./;)/(./i./2./3./;) x (4) 

JM 



,J+i3+i4-M-m3-m4^ imij3 _ ms J j _ M)(j 2 m 2 j 4 - m 4 I JM). 
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In the above equation, (jx\ \Yj\ [73) is the reduced matrix element of the spherical har- 
monics Y Jfl , J = \/2J + 1, and /(jij2j3i4) is the radial integral: 

roo 

I{hhhh) = J o M^Uj^Uj^Uj^Uj^r)—, (5) 

where u(r) is the radial part of the HF single-particle wavefunction. As it is shown in 
|13|14j the radial integrals can be calculated accurately by choosing a large enough cutoff 
radius R and using a iV-point integration Gauss formula with abscissas and weights r k ,w k - 

I{jihhh) - Nq 1 ^ u h {r k )u j2 {r k )u h {r k )u j4 {r k ) (6) 

Z k=l r k 

So we employ the hamiltonian including an average nuclear HF field, pairing interac- 
tions, the isoscalar and isovector particle-hole residual forces in the finite rank separable 
form |14j . This hamiltonian has the same form as the QPM hamiltonian with N separable 
terms |4ll8j . but in contrast to the QPM all parameters of this hamiltonian are expressed 
through parameters of the Skyrme forces. 

In what follows we work in the quasiparticle representation defined by the canonical 
Bogoliubov transformation: 

4» = «;<*t m + (-iy- m v jaj _ m . (7) 

The single-particle states are specified by the quantum numbers (jm) The quasiparticle 
energies, the chemical potentials, the energy gap and the coefficients u,v of the Bogoliubov 
transformations (JTj) are determined from the BCS equations. 
We introduce the phonon creation operators 

= \Y. ( X jj' A+ Uf> V) - (-1) A -"1# A(jj'; X - p)) (8) 
jj' 

where 

^ + (j/;V) = Yijmj'm I X/j)aj m at ,. (9) 
/ J 

mm 

The index A denotes total angular momentum and [i is its z-projection in the laboratory 
system. One assumes that the quasiparticle RPA (QRPA) ground state is the phonon 
vacuum | 0), i.e. Qx^i | 0) =0. We define the excited states for this approximation by 

Making use of the linearized equation-of-motion approach pQ one can derive the QRPA 
equations 01]: 

(10) 

In QRPA problems there appear two types of interaction matrix elements, the ^(jlj' )(j 2 j') 

matrix related to forward-going graphs and the &{jlj>)(j 2 j' 2 ) matrix related to backward- 
going graphs [Til . Solutions of this set of linear equations yield the eigen-energies and the 
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amplitudes X, Y of the excited states. A dimension of the matrixes A, B is a space size 
of the two-quasiparticle configurations. 

Using the finite rank approximation we need to invert a matrix having a dimension 
AN x AN independently of the configuration space size. One can find a prescription how 
to solve the system (fTU|) within our approach in |13ll4j . The QRPA equations in the 
QPM |4|18j have the same form as the equations derived within our approach[13 14, but 
the single-particle spectrum and parameters of the p-h residual interaction are calculated 
making use of the Skyrme forces. 

In this work we use the standard parametrization SIII |TjJ| of the Skyrme force. Spherical 
symmetry is assumed for the HF ground states. It is well known |10|llj that the constant 
gap approximation leads to an overestimating of occupation probabilities for subshells 
that are far from the Fermi level and it is necessary to introduce a cut-off in the single- 
particle space. Above this cut-off subshells don't participate in the pairing effect. In 
our calculations we choose the BCS subspace to include all subshells lying below 5 MeV. 
The pairing constants are fixed to reproduce the odd-even mass difference of neighboring 
nuclei. In order to perform RPA calculations, the single-particle continuum is discretized 
[To] by diagonalizing the HF hamiltonian on a basis of twelve harmonic oscillator shells 
and cutting off the single-particle spectra at the energy of 190 MeV. This is sufficient to 
exhaust practically all the energy-weighted sum rule. 

Our investigations enable us to conclude that N=A5 is enough for multipolarities 
A < 3 in nuclei with A < 208. Increasing N, for example, up to N=60 in 208 Pb does not 
change results for energies and transition probabilities practically. Our calculations show 
that, for the normal parity states one can neglect the spin-multipole interactions as a rule 
and this reduces by a factor 2 the total matrix dimension. For heavy nuclei our approach 
gives a large gain in comparison with an exact diagonalization [T4*] . 

To take into account the mixing of the configurations in the simplest case one can write 
the wave functions of excited states as: 



i AiiiA2«2 

with the normalization condition: 

(*„(jm)\mjm)) = J:rKJv) + 2 E (Px£(J")f = l (12) 

i AiuA2i2 

The matrix element coupling one- and two-phonon configurations is: 

(Q Jt \H\ [Qt in QU]j) = Uj${Ji) (13) 

U\lll{Ji) is some combination of the geometrical factors and the QRPA phonon ampli- 
tudes gED| ■ 

The energies of the states ^(A/i) are solutions of the following equation jl]: 



F(E U ) = det 



2 A x ii,Aa»2 UXlil + ~ Ev 



(14) 
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The rank of the determinant equals the number of the one-phonon configurations in- 
cluded in the first term of the wave function ty u (\fi). 

It is worth to point out that after solving the RPA problem with a separable interaction, 
to take into account the coupling with two-phonon configurations demands to diagonalize 
a matrix having a size that does not exceed 40 even for the giant resonance calculations 
in heavy nuclei whereas one would need to diagonalize a matrix with a dimension of the 
order of a few thousands at least for a non-separable case. 

3. RESULTS OF CALCULATIONS 

As an example we consider the 2^, 3f state energies and transition probabilities B(EX) 
in 112 Sn. The experimental data |21|22j and the results of our calculations within the 
QRPA (the second line) and with taking into account the two-phonon terms (the third 
line) are shown in Table 1. In our calculations the two-phonon terms of the wave func- 
tion include phonons with multipolarities A = 2,3,4,5. One can see that there is 
some overestimate of the energies and transition probabilities for the QRPA calculations. 
The inclusion of the two-phonon configurations results in a reduction of the enrgies and 
transition probabilities for the 2^,3^ states in 112 Sn. Generally there is a reasonable 
agreement between theory and experiment. The study of an influence of a choice for the 
Skyrme forces parameters on properties of the low-lying states within our approach and 
calculations for other nuclei are in progress now. 

Table 1 

Energies and B(EA)-values for up-transitions to the first 2 + , 3~ states in 112 Sn. 



State 




2t 




3f 




Energy 
(MeV) 


B(E2) 
(e 2 fm 4 ) 


Energy 
(MeV) 


B(E3) 
(e 2 fm 6 ) 


EXP. 


1.26 


2400±140 


2.36 


87000 ±12000 


QRPA 


1.49 


2600 


2.73 


97000 


2PH 


0.90 


2200 


1.90 


72000 



4. CONCLUSION 

A finite rank separable approximation for the QRPA calculations with Skyrme interac- 
tions that was proposed in our previous work is extended to take into account the coupling 
between one- and two-phonon terms in the wave functions of excited states. 

It is shown that the suggested approach enables one to reduce remarkably the dimen- 
sions of the matrices that must be inverted to perform structure calculations in very large 
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configuration spaces. 

As an illustration of the method we have calculated the energies and transition proba- 
bilities of the 2f and 3j~ states in 112 Sn. 

They are in a reasonable agreement with experimental data. A systematical study of 
an influence of the two-phonon terms on properties of the low lying states is in progress 
now. 
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